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I. INTRODUCTION 



Entanglement of a quantum system not only has important consequences as a fundamental physical property, it has 
also been identified as a resource in quantum information process (QIP), quantum computing and quantum cryptogra- 
phy [l|. Entanglement is fragile when the system of interest is not isolated from the influence of its environments. The 
dynamical aspect of entangled state in the context of quantum open system has been investigated in many different 
scenarios such as disentanglement of qubit systems , continuous variable systems 0, [|[ , entanglement delayed 
creation and revival and non-Markovian entanglement evolution [l3|, [TiJ , to name a few. 

As is well known, environmental noises that cause quantum systems to decohere also lead to the loss of entanglement 
p^| - [l7j . In its simplest form the system plus environment can be examined in different configurations of system- 
environment models including (i) harmonic oscillators coupled to a bath of bosons, i.e. Quantum Brownian motion 
model [HI, HH; (ii) Two-level atoms or spin-1/2 particles (Qubits) interacting with a bath of harmonic oscillators, 
i.e. spin-boson model 0, IH| and (iii) Qubits (or harmonic oscillators) coupled to a bath of spins-1/2, i.e. spin bath 
model |2j2j-[33j. It is shown that open system dynamics can be significantly modified by environmental configurations. 

In this paper, we discuss entanglement dephasing dynamics for spin environment models. Our strategy is to directly 
solve the models without any approximations that ma y le ad to the conversion of bath spins into the effective Holstein- 
Primakoff bosons or Schwinger bosons (e.g., see j23ll25j). The purpose of this paper is to find important features 
induced from the spin environment that are not seen in the boson case. We not only consider different couplings 
between the system and its environments, but also take into account the effect of the various initial states of the 
environment on the entanglement dynamics of the central qubit system. Here we take our central qubit system as 
two uncoupled spin-1/2 particles (qubits), which are interacting with a bath of spins-1/2 particles through z-x type 
coupling. We have presented several analytical results on the dynamics of the two central qubits in physically relevant 
situations including individual coupling model and collective coupling model. 

This paper is organized as follows. In Sec|TT]we present the Hamiltonians of two-qubit dephasing models with local 
and global spin environments. The analytical solutions to the both models are provided. The entanglement dynamics 
of the central qubits with various initial states are analytically investigated in Sec lIIII for the local coupling, and in 
SeclIVI for the global coupling, respectively. Specifically, for each model we consider three different types of initial 
states of the environment: pure, mixed and entangled. In addition, we have extended our analytical treatments to 
numerical simulations for more generic situations in SecJVl it is shown that the concurrence dynamics of the central 
qubits is sensitively dependent on the total numbers of environment spins and bath initial states. We conclude in 
SecEIl 



II. MODELS AND SOLUTIONS 



We consider two well-known cases for system and environment coupling: Case (a), the two system qubits are 
coupled to a common environment; Case (b), the two system qubits are interacting only with their local environments 
and there is no crosstalk between the two environments. For simplicity, the interaction between spins is neglected. 
Precisely, the total Hamiltonian for these two cases are given by: 

Htot = H sys + H cnv + -Hint) (1) 
H sys = u A cr^ + uj B af , (2) 

where for the global environment case we have H CDV — Ylk u}k<J z> ^ int = + a z)^2k Ck<J x- ^ n * ne case °f l° ca l 
environment we have H cnv = J^ fc u>k&z + J ^2 l 0Jia l z , H- lrA — ^2 k Ck<r x + erf ci<r l x . oja and ojb are the transition 
frequencies for the two qubits, respectively. For simplicity, but without loss of generality, we assume that uja = = 
lj s . Here, a z and a x are the Pauli matrices, and Uk is the frequency for fcth spin in the environment and Ck are 
dimensionless real- value coupling constants between the feth spin and the central qubits. As is shown in Fig. []Ja), ^ 
and I represent the two local baths for the central qubits respectively. In both cases, the dephasing coupling ensures 
that quantum coherence is damped, but the energy of central qubits is perfectly preserved. 

For both cases, the dynamics for the central qubit systems can be solved analytically under the assumption that 
initially the system and environment are in a separable state, i.e. the density matrix for the system plus environment 
at t — is of the following form: 

p(0) = Psys (0)®p B (0). (3) 

The environment spins may take a simple form ps(0) = rife/°s(0)- Otherwise, the environment may contain some 
entangled spin blocks, in this case, the environment state may be written as /5b(0) = Y[j Ps(0): where p J B stands for 
the density matrix for each block (More precise descriptions are given in Subsection IIII Dj) . 
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(a) 

Local Environment 

(b) Global Environment 




FIG. 1. (a) Two qubits are interacting with their local environments. Two central qubits are prepared in an initially entangled 
state but have no interaction between each other, (b) Two qubits are interacting with one global environment. 



For our dephasing model, the time-dependent density matrix for the total system may be written as: 

pit) = e ~ iHMt p(0)e iHtatt = e-^^tJJ e -4 **' p(0) TJ e ilfct ]e^*, (4) 

k k 

where -ff sys is the system Hamiltonian and Xk is the sum of the bath Hamiltonian and interaction Hamiltonian of the 

fcth environment spin. For example, in case (b), Xk = oJkO'z + {°~z + <T ?) c fc cr x- The elements of the reduced density 
matrix (RDM) for the central qubit system can be obtained by tracing over the bath spin degree of freedom: 

p%°(t) = (m\Tr k p(t)\n) =6^-^(0)11^4 = ^^-^'fifO)/™!*), (5) 

k 

f mn (t) = n Tr ^™« ^l[Tv k [e-^ E ^p k (0)e^ E ^}. (6) 

k k 

\m), \n) and E' m , E' n are the eigenvectors and eigenvalues of the system Hamiltonian H sys respectively. Note that the 
notations E m ,E n denote the eigenvalues of the coupling operator a^ + a B , thus E' mn = u> s E m , n . It should be noticed 
that the dephasing coefficient f mn (t) is a time dependent function satisfying |/ mn (t)| < 1 and f nn = 1. The specific 
form of f Trinity depends on the interaction Hamiltonian and the initial state of environment, which will be discussed 
in later sections. 



III. SOLUTION OF LOCAL COUPLING MODEL WITH DIFFERENT INITIAL STATES 



A. General Solutions 



When each central qubit is interacting with its own environment, as shown in Fig. [Ha), the RDM of the two central 
qubits can be obtained by using Kraus operator technique [5|,l3J]- In principle, we can always write the time-dependent 
RDM in such a way: 

p(t)=Y,K iP (0)Kl (7) 

i 

where Ki satisfy the condition ^\ K\Ki = I. It is convenient to solve the model in this way since we may just focus 
on the subsystem involving only one qubit plus its local environment. The RDM of such a subsystem is thus a 2 x 2 
matrix, which can be derived directly from Eq. ([5]): 

Obviously for the dephasing model, the diagonal elements are time independent, i.e. p n n(t) — p n n(0), so we only 
need to consider the off-diagonal ones, for which E m E n = — 1. Since pi2 = p\ x for a 2 x 2 density matrix, we define 
f(t) = fi2(t)- The Kraus operators of one-qubit subsystem are given by: 



V - I 1 ) K A ( ° 



(9) 
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where f(t) is the dephasing coefficient defined in Eq. ([5]). Then the RDM for two qubit system can be directly obtained 
by: 



p AB (t) = £ K AB p A *(0)K 
»,i=i,2 



AB-f 
ij ) 



where the composite Kraus operators Kij = K A ® Kf . are given by 



K 21 
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Then for an arbitrary initial state, the time dependent RDM is simply given by: 

/ Pn e- i2uAt f(t)p 12 e~ a " At f{t)p 13 e~^ At f{tfp 1A \ 

, A _ e^ t f(t)*p 21 p 22 f(t) 2 p 23 e- i ^ t f(t)*p 2i 

m ~ e^f(ty P31 f(ty 2 P32 P33 e^ At f{ty P3A 

\e^ At f{t)* 2 Pil e i2 ^*/(i)p 4 2 e a " At f{t)p i3 p 4i J 



(11) 



Note that the solution (fTTj) is obtained without any approximations. We will discuss three special cases where the 
function f(t) can be explicitly evaluated. 



B. Mixed State Environment 

As our first example, we assume that the environment spins are in the following mixed state: 

o k (0)-( N + ° 
P W - { o N k 



(12) 



where N k ,N k are the probability of which the fcth spin is in up or down state (N!L + N't — 1). Here, we may choose 



/(e^ fc + e P u ' k ) 1 where (3 = 1/T and T is the temperature. By using the following identity 



i sin \J u 2 + v 2 
VvF+v 2 



(ua z + va x ) 



and Eq. 



it is easy to show that: 
Tr fe [0 = cos 2 (p fe i) + 



sin (Pkt) , 2 . v v 2-, fAT k, AT ks , 2/-, P P ^Sinjpfet) 
2 (w fc + -E m -E„c fe )(A' + + A'_) = 1 -c k (l -E m E n ) 5 , 



(13) 



where p k = \/^l+ cj- 
From (fl~3|) . we obtain: 



Pi A 2 W = Pi 4 2 (o) e - j(B ™^" ) 7W, 



/w=n[i-2ci 



'Ml 



(14) 
(15) 



Eq. (|15[) is the exact solution in the condition of finite number of environment spins. It is easy to see that each 
environment spin contributes an oscillating function, which varies from —1 to 1, affecting the coherence property of 
the central qubit system. When the environment involves more and more spins, the overall effect of all spins results in 
a typically fast decaying function f(t). It is evident from Eq. (fT3")) that the RDM is insensitive to the distribution of 



5 



. However, it is worth noting that this is not a generic feature for the spin environment. Actually, as shown 
in subsection IIVB1 the coherence function f(t) can explicitly depend on the distributions N^,N^ . 

Interestingly, if more information about the bath spectrum density becomes available, we may get a rather compact 
form for the function /(£). For example, for a bath with infinite numbers of spins and the Ohmic spectrum density, 

J M = V-^e'^, (16) 

with the weak coupling approximation Cfc <C 1, we can get: 

/(t) = (l + 4u, c 2 i 2 )-H (17) 

The power law decay of quantum coherence clearly deviates from the exponential decay commonly seen in the case 
of Markov bosonic bath [16[ . 



C. Pure State Environment 



Another interesting case is that the bath spins are initially prepared in a pure state: 

|*(0)) = cosa fc | t) +sina fc | |). 

The density matrix may be written as, 

/(0) = |*(0))<*(0)| = ( COS "\f ak ) ■ (18) 

V cos ckfc sin ctk sin ak I 

The RDM for the central qubits at time t is also determined by the same terms as in Eq. @, thus we have 

2 

Tr fe [F„ fc m ] = 1 - % sin 2 (p fc <)(l - E m E n ) + i^- sm{p k t) sm2a k (E n - E m ). (19) 

Eq. (j!9[) is different from Eq. (|13[) by an imaginary factor in the order of c k /p k . Finally, the time dependent dephasing 
coefficient is given by 

2 

/(*) = TTl 1 - 2%sin 2 (p fc i) sin(p fc i)sin2Q fc ]. (20) 

k Pk Pk 

We note that the above expression does not permit fit) to have the same compact and useful form as Eq. (IT71) for 
the continuous spectrum density considered before. 



D. Environment with Entangled Blocks 



We also consider a situation where the environment spins may be decomposed into a set of entangled blocks. This 
clearly is an interesting configuration that may impose some more strong constraints on the central qubit entanglement 
evolution. More specifically, we begin with an environment consisting of TV x M spins in total. We may divide the 
total N X M spins into M blocks, each contains exactly N spins, which is prepared in an iV-GHZ state [35|. In this 
case, the matrix elements of the RDM are given by Eq. (|5J with a modified dephasing function f mn {t): 

UN N M 

a) - n t ^ i n e-^p^m n ^ = n ^K^m 

j=l k=l k=l j=l 

where pg J (0) is the density matrix for entangled bath spins confined in the jth block. And we have 

N 9 N 

rp ipNJm TJn c k ■ 2 r, p», / E n - E m N 1 - cos 2p k t sm2p k t 

TTF mn (t) = I 1(1 2 Sm ~ E m E n )) + ) UJ k Ck + I Cfe) 

£Ji Pk 2 11 p\ pk 

N 

E n -E m . NTJ l-cos2p fc t ,sm2p k t 

+ 5 II 2 U k C k +l Cfe. (21 

2 £Ji Pt Pk 
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A particular expression for the coherence function f(t) may be obtained when each block only contains two entangled 
spins, which are prepared as 3^' = ^=(| H) + | tt)) and Cji = Cj2 = Cj,pji = pj2 = Pj, then the off-diagonal terms 
in the RDM are: 

P 3 m ^n(t) = e-^ E -- E ^p s m ^ n (0)f(t), (22) 

(23) 



M 



f(t) = Y[(l-2^ S m 2 2 Pj t). 

Again, with the Ohmic spectrum density assumption and c k « l,the dephasing coefficient f(t) in the long time 
limit is simply given by: 



/(i) = (l + 16 Wc 2 i 2 )-^«( 4Wc <)-". 



(24) 



Immediately, we can see from the above calculations that the block-entangled environment may give rise to more 
severe decoherence effect on the central qubits than the mixed initial state (e.g., see Eq. (|T7| ). It is interesting to 
consider another simple case: Cji — —Cj2 = Cj,Pji = Pj2 = Pj, with which we get 



M 

/w=ri( i 



2 2 

5^/sin 4 ^t . 
Pi 



(25) 



Using the same assumption and approximation above, we get a simple function in the long time limit: 

f(t) = (1 + 16oj 2 t 2 )^(l + 4w 2 t 2 )- 2 " « (4w 2 t 3 )-". (26) 

IV. SOLUTION OF GLOBAL COUPLING MODEL WITH DIFFERENT INITIAL STATES 

A. General Solution 

In this subsection, we consider the case that both qubits are coupled to the same environment, illustrated in 
Fig. [T](b). Again, this model can be solved exactly by using the same method as we used in the previous section, 
but with different conventions E m ,E n 6 {2,0,0,-2}. It can be seen below that new features can arise from the 
backreaction on the central qubits induced by the common bath. Thus, the resultant RDM is different and more 
complex, which is given by: 



/ Pn 

p 2 ie i2 "° t g(ty 



Pi2e- i2 "° t g{t) Pi^^git) p 14 e 



P22 
P32 



V PAie^fit)* p^e^git) 



P±3C 



P23 
P33 



P24e 
P34e' l2 ^ t g(t) 

P44 



(27) 



/ 



for simplicity. Considering the the initial states of system qubits A, B are in Bell states. When ^ AB = -J=(| ft)±l U)), 
the result RDM will be similar as in Eq. flTTJl with f(t) 2 -t f(t). Yet for the other two states ^ AB = ^(| fl) ± | It)), 
it can be shown that entanglement is perfectly preserved. 



B. Mixed State Environment 



For the environment spins that are initially in a mixed state, we have 
f(t)=H[l-8^ S m 2 (q k t)}, 

k q k 

i \ TT f / \ / \ sin(<7fci) sin(cjfct) 
9{t) = M cos(g fe i) cos(u k t) + uj k h iA k 

where q k = ^jf+Ic 2 , A k = N i +N ^ . 



(28) 



cos(q k t) sin(w fe t) - — sin(q k t) cos(uj k t) 
qk 
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C. Pure State Environment 



If the bath is initially in a pure state (see Eq. (fT8l) . by assuming a k — 7r/4), we may get an explicit form of the 
density matrix Eq. (|27p with 



2 

/(*) = TTi 1 - 8 % s ™ 2 (qkt) + i2— sm[2q k t}}, (29) 
k q k qk 

g(t) = TT{cos((7 fe i) cos(u} k t) H sin(<j fc i) sin(w fe i) - i— sin(q k t) cos(wfei)}. 

" 9fc <7fc 

This solution will be used for entanglement dynamics to be discussed later. 



D. Environment with Entangled Blocks 

For simplicity, we consider a case in which the environment has M blocks, each block consists of two entangled 
spins (i.e., N = 2). For this structured environment, the time-dependent coefficients in Eq. ([27)) takes the following 
form: 

M 

f = \[ A - 84 1 5| >1 )(1 - 8c£ 2 S£a) + Wc^Cj^i^S^Sli ~ C JtL C^S jA Sj fi )], (30) 
o 

M 

9 — JJ[(Cj,ioCj,i + ^j,io<S'i,iw| i i)((7 7 - i 2oC^,2 + SjtwSj^Wjfi) 
j 

— ujj^ujj^iCjsoSj^ — Sj t ioCj t i) (Cj^oSj^ — Sj t 2oCj,2) — ^Cj^Cj^Sj^Sj^ cos^uij^ +^,2)]- (31) 
with the notation 

Cj ti = coa(qj li t),Sj t i = sm q ^ %t ; c jiQ = cos{ojj fi t),Sj,iQ = s\n(uj hl t),q .i = J~tjf ti + 4c^, i = 1,2. 
where w^j and Cj^ are the frequency and coupling strength for the «th spin in the jth block. 

V. ENTANGLEMENT DYNAMICS FOR SPIN ENVIRONMENTS 

In this paper we use Wootters' concurrence to measure the entanglement of the central qubits [36| . which is defined 

as 

C(p) = max{0, VXi - n/a 2 - VX 3 - V^ 4 }, (32) 
where Ai are the eigenvalues in decreasing order of the matrix 

C = P((7y ® <7y)p*(<7y ® (Jj,). 

For our model all the off-diagonal elements in the RDM are obtained by their initial value times a dephasing coefficient, 
thus if any element is zero at t = 0, then it will remain zero for any t > 0. Consequently, when Bell states are chosen 
as the initial state for the central qubit system, its RDM will preserve an X-form [34|]. In this case, the concurrence 
may be simplified as: 



C(p) = 2max{0, |p 14 (t)| - VrfW), \paa(t)\ - V 'pu(t) Pii (t)} = 2max{0, \f{t)\}. (33) 

It is clear that the local environments do not generate entanglement between central qubits, so our focus will be 
on entanglement decay. For the global environment, however, it is possible to achieve entanglement generation when 
initially the central qubits are in a separable state. Based on our analytical solutions to both models, we will be 
able to systematically study concurrence dynamics C of the central qubits in the following two subsection. In all the 
numerical results presented below, we will make the following conventions: the frequency of system qubits is chosen as 
lo s = 1; for the environment spins, their frequencies and coupling constants are randomly distributed in the intervals 
uj k E (1,2), c fe € (0.1,0.2), respectively. 
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FIG. 2. Two qubits interacting locally with two 6-spin environments. The blue solid line is for environment initially in mixed 
state, red doted line for pure state, and yellow dashed line for block-entangled environment. 




FIG. 3. Two qubits interacting locally with two 60-spin environments. 
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FIG. 4. Two qubits interacting locally with two 400-spin environments. 



A. Entanglement Decay in Local Environment Model 



Now we consider the central qubit system to be prepared in the Bell state ^ ab — (| tt) + I ID)- Figs. [2J [3] and 
2] show the concurrence dynamics with 6, 60 and 400 environment spins. Here and next section, the blue solid line is 
always used for the situation that the environment is in the mixed state (for local environment Afc is not important, 
but we have set = 0.9 in next subsection for global environment.), the red dotted line is for environment in the 
pure state and the yellow dashed one is for block-entangled environment. 

If the number of the environment spins is not too large, as illustrated in Figs. [5] and the central qubit entanglement 
is a quasi-periodic function of i, representing a strong non-Markovian feature of environment since the information 
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FIG. 5. Two qubits interacting with a 6-spin global environment. The blue solid line is for environment initially in mixed state, 
red dotted line is for pure state, yellow dashed line is for pair-wise entangled bath. 
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FIG. 6. Two qubits interacting with a 60-spin global environment. 

transferred to the environment from the system will quickly go back to itself. It should be noted that block-entangled 
environment could play a positive role in preserving the system concurrence compared with the cases of mixed or pure 
initial states. In general, the mixed environment state may cause more detrimental effect to entanglement of central 
qubits. 

Typically, a larger environment shows more dramatic differences between the three initial states. For the 400- 
spin bath shown in Fig. 21 we see that the initial entanglement decays to zero much quicker than the case of small 
environments. This is consistent with our general observations for the bosonic case. In all the three cases, the increased 
degrees of freedom of the environment will cause effectively an irreversible information flow between the system and 
its environment. It is noticed that concurrence revival is found in the case of both the block-entangled environment 
and pure state, but not in the mixed state case. For our example, as shown in Figs. [5] and [31 entanglement decay and 
revival occur approximately at the same pace for mixed and pure initial states. But for the block-entangled bath, the 
structured environment may speed up both decay and revival processes. 



B. Entanglement Generation in Global Environment Model 

When the system qubits are coupled to a global environment, it is possible to achieve entanglement generation 
between the two central qubits. It is clear that such generation of entanglement is partly due to the back-reaction 
on the central qubits through system-environment couplings. For example, let us consider a separable initial state 
of the central qubits A, B represented by ^ ab — \{\ \)a + I i) a) ® (| t) b + | \)b)- In what follows we show that 
entanglement generation is sensitively dependent on the environment initial states. 

We begin with the case of small environment with only 6 spins, as shown in Fig. [5] We see again a quasi-periodic 
evolution as expected for a small environment. A more interesting phenomenon is the generation of entanglement 
by the mixed state bath. We have found that the mixed initial state is very efficient in producing high degree of 
entanglement. Moreover, it shows that concurrence proceeds in time in a predictable manner for a long period of 
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FIG. 7. Two qubits interacting with a 200-spin global environment. 
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FIG. 8. Two qubits interacting with a 800-spin global environment. 

time. Although the pure state environment can also generate entanglement, but the degree is much lower than the 
mixed state. For the case of the block-entangled environment (it is assumed each pair is in a GHZ state), surprisingly 
we find that, for possible configurations, there is no generation of entanglement. 

If the environment contains a large number of spins, the pure state environment may generate a high degree of 
entanglement. On the contrary, for the mixed state environment, less entanglement is generated. Fig. [6] shows the 
entanglement dynamics for two central qubits coupled to a global bath with 60 spins. The entanglement generated may 
only survive some time before decaying to zero, and stay as zero for a long period of time (at least no revival is found 
in our simulations). However, the entanglement generated by the mixed state environment shows a regular creation 
and decay pattern that is not seen for the pure initial state. Similar oscillating patterns are also observed when the 
number of spins is increased to 200. In Fig. [3 two bumps are found at early stage, and no visible entanglement is 
detected at a late stage. This phenomenon occurs when the environment is larger than the size of system, but not 
large enough to trap all the information leaked into the environment from the system. 

Finally, we plot Fig. [8] for a very large environment consisting of 800 spins. In many ways, the model under 
this condition may be well represented by a Markov system. Intuitively, for such a large environment, the system 
information once lost into the environment won't come back to itself again. For our example, entanglement generation 
is not observed for the mixed state bath. But for the pure state bath, we see 0.4 entanglement generation measured 
by concurrence in an early stage due to the initial relaxation of the environment into its steady state. Such an initial 
generation of entanglement is also observed in the Holstein-Primakoff boson case [23|, [25| . 

VI. CONCLUSION 

We studied dynamics of entanglement of the open system under the z-x interaction with different spin environments. 
The all models presented here are solved analytically. The power law decay of entanglement is obtained when the 
environment is large and the spectrum density is of an ohmic type. 

For the local environment model, we found that, in the case of the small environment, the central qubits can 
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typically maintain a high level entanglement for a long period of time. Interestingly, we have shown that the block- 
entangled environment provides stronger protection of the entanglement. In addition, the entangled spins confined in 
the local environments can accelerate the decay and revival of the entanglement between central qubits. 

For the global environment model, we found that, in the case of small environments (see Fig. [5]), the concurrence 
evolves in a similar pattern when driven by a pure state or mixed state baths. Even for a large environment (see, Fig. [6] 
and Fig. [7]), the entanglement generated is still seen to decay in a non-monotonic way, a feature that is commonly 
ascribed to the non-Markovian dynamics. For a very large environment (see, Fig. [5]), similar to Markov limit, the 
mixed state environment considered here does generate entanglement. However, the pure state bath entanglement 
generation at early stage is possible. As for block-entangled bath with GHZ state as the initial state, the generation 
of entanglement is always suppressed in the examples considered here; a more detailed discussion of this interesting 
phenomenon will be discussed in another paper. 
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